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Context. X-ray observations of evolved supernova remnants (e.g. the Cygnus loop and the Vela SNRs) reveal emission originating from the 
interaction of shock waves with small interstellar gas clouds. 
, Aims. We study and discuss the time-dependent X-ray emission predicted by hydrodynamic modeling of the interaction of a SNR shock wave 
with an interstellar gas cloud. The scope includes: 1) to study the correspondence between modeled and X-ray emitting structures, 2) to explore 
two different physical regimes in which either thermal conduction or radiative cooling plays a dominant role, and 3) to investigate the effects 
of the physical processes at work on the emission of the shocked cloud in the two different regimes. 
Q ' Methods. We use a detailed hydrodynamic model, including thermal conduction and radiation, and explore two cases characterized by different 
?—j | Mach numbers of the primary shock: A4 = 30 (post-shock temperature r ps h ~ 1.7 MK) in which the cloud dynamics is dominated by radiative 
C/2 . cooling and M = 50 (7" psh as 4.7 MK) dominated by thermal conduction. From the simulations, we synthesize the expected X-ray emission, 
using available spectral codes. 

^> Results. The morphology of the X-ray emitting structures is significantly different from that of the flow structures originating from the 
• '■H , shock-cloud interaction. The hydrodynamic instabilities are never clearly visible in the X-ray band. Shocked clouds are preferentially visible 
' during the early phases of their evolution. Thermal conduction and radiative cooling lead to two different phases of the shocked cloud: a cold 
cooling dominated core emitting at low energies and a hot thermally conducting corona emitting in the X-ray band. The thermal conduction 
makes the X-ray image of the cloud smaller, more diffuse, and shorter-lived than that observed when thermal conduction is neglected. 
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1. Introduction 

This paper is part of a series devoted to study the interaction 
of a shock wave of an evolved supernova remnant (SNR) with 
a small interstellar gas cloud (like the ones observed, for in- 
stance, in the Cygnus loop or in the Vela SNR) through detailed 
and extensive numerical modeling. The project aims at over- 
coming some of the limitations found in the previous analogous 
studies (for instance, taking into account simultaneously im- 
portant physical effects such as heat flux and radiative losses) 
and crucial for the accurate interpretation of the high resolu- 
tion multi-wavelength observations of middle-aged SNR shell 
available with the last-generation observatories. 
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In a previous paper dOrlando et alJl2005l hereafter Paper 
I), we have modeled in detail the shock-cloud interaction with 
hydrodynamic simulations including the effects of thermal con- 
duction and radiative losses from an optically thin plasma. We 
have investigated the interplay of the latter two processes on 
the cloud evolution and on the mass and energy exchange be- 
tween the cloud and the surrounding medium, by exploring two 
different physical regimes in which one of the two processes is 
dominant. We have found that, when the radiative losses are 
dominant, the shocked cloud fragments into cold, dense, and 
compact filaments surrounded by a hot corona gradually ab- 
lated by the thermal conduction; to the contrary, when the ther- 
mal conduction is dominant, the shocked cloud evaporates in 
a few dynamical time-scales. In both cases we have found that 
the thermal conduction is very effective in suppressing the hy- 
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drodynamic instabilities that would develop at the cloud bound- 
aries. 

In this paper, we study and discuss the time-dependent 
X-ray emission predicted by the modeling mentioned above. 
Several authors have investigated the emission from SNRs in 
different spectral bands, adopting global SNR models includ- 
ing the effects of radiative cooling, evapor ating clouds, large 



Table 1. Summary of the initial physical parameters character- 
izing the unperturbed ambient medium and the spherical cloud. 
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On the other hand, the X-ray emission originating from model 
shocked interstellar clouds has not been investigated yet in de- 
tail, despite it is potentially important for the energy budget of 
the shocked ISM and for the interpretation of the observations. 

Here we synthesize from the numerical simulations de- 
scribed in Paper I the X-ray emission expected from the shock- 
cloud interaction. Our scope includes: 1) to link modeled to 
X-ray emitting structures, 2) to investigate the emission of the 
shocked cloud in two different physical regimes in which either 
thermal conduction or radiative cooling is dominant, and 3) to 
investigate the effects of thermal conduction and radiation on 
the emission of the shocked cloud. 

The paper is structured as follows: in Sect. |2] we briefly 
summarize the model of the shock-cloud collision and outline 
the method to synthesize the X-ray emission from the numeri- 
cal simulations; in Sect.[3]we discuss the results; and finally in 
Sect.|4]we draw our conclusions. 

2. The modeling 

2.1. Hydrodynamic simulations 

In this section, we summarize the model of the shock-cloud 
collision. We refer the reader to Paper I for more details. 

The model describes the impact of a planar shock front 
onto an isolated gas cloud. The cloud before the impact is as- 
sumed to be spherical with radius r c \ = 1 pc, small compared 
to the curvature radius of the SNR shock 1 , internally isother- 
mal, and in pressure equilibrium with the surrounding medium. 
The unperturbed ambient medium is assumed to be isothermal 
(with temperature Ti sm = 10 4 K, corresponding to an isother- 
mal sound speed c s = 11.5 km s~') and homogeneous (with 
hydrogen number density «j sm = 0.1 cm 4 ). The total mass of 
the cloud is ~ 0.13 M sun . Table^summarizes the initial physi- 
cal parameters characterizing the unperturbed ambient medium 
and the spherical cloud. The shock propagates with a velocity 
w = A1cj sm in the ambient medium, where Al is the shock 
Mach number, and c; sm is the sound speed in the interstel- 
lar medium; the post-shock conditions of the ambient medium 
well before the impact onto the c loud are given by the strong 
shock limit (Zel'dovich & Raizer 1966, see also Paper I). The 
fluid is assumed to be fully ionized, and is regarded as a perfect 
gas (with a ratio of specific heats y = 5/3). 

The plasma evolution is derived by solving the time- 
dependent fluid equations of mass, momentum, and en- 
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ergy conservation. We take into account the thermal con- 
duction (according to the formulation of Spitzer 1962), m- 
clud ing the free-streaming limit (saturation) on the heat flux 
JCowie & McKeel 19771 [Giuliani 1 9841 iBorkowski et alJI 19891 
Fad evev et al.ll2.Q02l and references therei n), and the radiative 
losses from an optically thin plasma (e.g. iRavmond & Smith! 
ll977[lMewe et al.l ll985 and later upgrades). Continuity equa- 
tion of a tracer of the original cloud material is solved in addi- 
tion to our set of hydrodynamic equations. The calculations are 
performed using the flash code dFrvxell et alJ l2000) with cus- 
tomized numerical modules that treat thermal conduction and 
optically thin radiative losses (see Paper I). 

To study the X-ray emission expected during the shock- 
cloud collisions, simulated in Paper I (with an effective spa- 
tial resolution of 132 zones per cloud radius) that allow 
us to explore two different physical regimes in which either 
thermal conduction or radiative cooling plays a dominant role. 
The set of simulations includes: two models neglecting ther- 
mal conduction and radiation and considering the At = 30 
and M — 50 shock cases in a 3-D cartesian coordinate sys- 
tem (x,y,z) (runs HYm30cl0 and HYm50cl0, respectively 2 ); 
two models with thermal conduction and radiation, considering 
the At — 30 and Al — 50 shock cases, and in a 2-D cylindri- 
cal coordinate system (r,z) (runs RCm30cl0 and RCm50cl0, 
respectively). The models neglecting both thermal conduction 
and radiation have been computed in a 3-D Cartesian coor- 
dinate system (x, y, z), in order to describe accurately the hy- 
drodynamic instabil ities developing a t the boundaries of the 
shocked cloud (e.g. IXu & St one 1995; see also Paper I). As 
we have demonstrated in Paper I, however, the heat conduction 
rapidly damps the hydrodynamic instabilities and, in this case, 
the essential evolutionary features of the system can be ade- 
quately captured in a model using 2-D cylindrical coordinate 
system (r,z). In all the cases considered, the cloud is initially 
10 times denser than the surrounding medium (hydrogen num- 
ber density of the cloud, n c \ = 1 cm 4 , see Tab. Table [2] 
summarizes the physical parameters characterizing the simula- 
tions, namely the shock Mach number, At, the density contrast 
between the cloud and the ambient medium, x — ici/«ism, the 
velocity of the SNR shock, w, the temperature and density of 
the post-shock ambient medium, T ps h and n ps h respectively, and 
the cloud crushing time, r cc , i.e. the characteristic tim e for the 
transmitted shock to cross the cloud JKlein et all 1994 see also 
Paper I). 



1 In the case of a small cloud, the SNR does not evolve significantly 
during the shock-cloud interaction, and the assumption of a planar 
shock is justified (see also K lein et alll994l) . 



2 Run HYm30cl0 has been derived from run HYm50cl0 through 
the scaling / — > tM, u — > u/M, T — > T/M 2 (where t is the time, 
u the gas velocity, and T the temperature), with distance, density, 
and pre-shock pressure left unchanged (the so-called Mach-scaling; 
iKlein et all 1994 see also Paper I). 
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Table 2. Parameters of the simulated shock-cloud interactions. 
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[km s~'] 


[10 6 K] 


[cm" 3 ] 
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HYm30clO' <; 


3-D cart. (x,y,z) 
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1.7 


0.4 


9.1 


no 


HYm50clO 


3-D cart, (x, y,z) 


50 


10 


574 


4.7 


0.4 


5.4 


no 


RCm30clO 


2-D cyl. (r,z) 


30 


10 


344 


1.7 


0.4 


9.1 


yes 


RCm50clO 


2-D cyl. (r,z) 


50 


10 


574 


4.7 


0.4 


5.4 


yes 



" Shock Mach number. ' Density of the post-shock ambient medium. 

* Density contrast cloud / ambient medium. ' Cloud crushing time (Klein et al. 1994). 

c Velocity of the SNR shock. s Run derived from HYm50cl0 through Mach 

d Temperature of the post-shock ambient medium. scaling (see Paper I). 



2.2. Synthesis of the X-ray emission 

From the model results we synthesize the X-ray emission of the 
shock-cloud system in different spectral bands of interest. The 
results of numerical simulations are the evolution of tempera- 
ture, density, and velocity of the plasma in the spatial domain. 
In the case of 2-D simulations, we reconstruct the 3-D spatial 
distribution of these physical quantities by rotating the 2-D slab 
around the symmetry z axis (r = 0). The emission measure in 
the j-th domain cell is emj = nj^Vj (where nL is the hydrogen 
number density in the cell, and V s is the cell volume). We de- 
rive distributions of emission measure vs. temperature, EM(T), 
in selected regions by binning the emission measure values in 
those regions into slots of temperature; the range of tempera- 
ture [4 < log T(K) < 7] is divided into 75 bins, all equal on a 
logarithmic scale. From the EM(T) distributions, we synthesize 
the X-ray spectrum, usin g the MEKAL spectral synthesis code 
Jlvlewe et alJl985tlKaastrall992land later upgrade s), assuming 
solar metal abundances ( Gre vesse & A nders 199JJ). 

To derive spatial maps of the X-ray emission from the 
shock-cloud system, we assume that the primary shock front 
propagates perpendicularly to the line-of-sight (in the follow- 
ing LoS) and that the depth along the LoS is 10 pc with inter- 
cloud conditions for the medium outside the numerical spatial 
domain. The X-ray spectra integrated along the LoS and on pix- 
els of size corresponding to the spatial resolution of the numer- 
ical simulations are then integrated in selected energy bands, 
obtaining the X-ray images of the shock-cloud system. 

3. Results 

In Paper I, we have studied and discussed the hydrodynamics 
of the shock-cloud interaction for the cases considered here. 
We found that the shocked cloud evolves in cold, dense, and 
compact cooling-dominated fragments surrounded by a hot di- 
luted thermally conducting corona, when the radiative losses 
are dominant (Ai = 30 shock case; see Figs. 7 and 8 in Paper 
I). In this case, the radiative cooling strongly modifies the struc- 
ture of the shock transmitted into the cloud, leading to a cold 
and dense gas phase. When the thermal conduction is the dom- 
inant process (Ai = 50 shock case), the shocked cloud evap- 
orates in a few dynamical time-scales, and a transition region 
from the inner part of the cloud to the ambient medium is gen- 
erated (see Figs. 4 and 5 in Paper I). 



3.1. Emission measure vs. temperature 

We use the cloud tracer mentioned in Sect.|2]to identify zones 
whose content is the original cloud material by more than 90%. 
From these zones, we then derive the EM(T) distribution of the 
cloud. Fig.[2shows the cloud EM(T) evolution for the Ai — 30 
(left panels) and Ai — 50 (right panels) cases, either without 
thermal conduction and radiation (hereafter HY models; up- 
per panels) or with both effects (hereafter RC models; lower 
panels). We show the EM(T) distributions sampled at steps of 
0.4 t cc since t - 0.2 t cc . 

Fig- [D shows that, in HY models, the EM(T) distribution 
of the cloud is steadily centered around the temperature of the 
shock transmitted into the cloud, T sc i » fiT^IX-* where (3 ~ 1.7 
(see Paper I), x = 10; we obtain T sd « 0.3 MK for Ai = 30 
and T sd * 0.8 MK for M = 50 (see dotted lines in Fig.[T}. The 
evolution of the EM(T) distribution for Ai = 30 a nd for Ai — 
50 ca ses is similar, according to the Mach-scaling JKlein et alJ 
1994, Paper I): it rapidly becomes quite broad, covering more 
than a decade in temperature around r sc i; then, at late stages, it 
gets narrower. 

The EM(T) distribution obtained from RC models sig- 
nificantly changes, depending on which process is dominant. 
When the radiative losses dominate (RCm30cl0; lower left 
panel in Fig. the EM(T) distribution below 1 MK evolves 
toward a steep power law (with negative index), drifting to the 
cold side due to the progressive cooling of the plasma. A small 
fraction of cloud material gradually thermalizes by conduction 
to the temperature of the surrounding medium, forming a small 
peak centered at T ps h- Thus, at variance from pure hydrodynam- 
ics, the plasma splits into two separate thermal components: a 
cold dense core (T < 0.1 MK, see Paper I) and a hot diluted 
corona (7" * r psh = 1.7 MK). 

In the conduction-dominated Mach 50 case (RCm50cl0; 
lower right panel in Fig.[0, the EM(T) distribution is initially 
broad and centered at the temperature of the shock transmitted 
into the cloud, r sc i; then its maximum gradually shifts to higher 
and higher temperatures up to r ps h ~ 4.7 MK getting more 
peaked, due to the thermalization of the cloud material to T p ^- 

3.2. X-ray emission 

Fig-Elshows the cloud X-ray light curves in the broad [0. 1-10] 
keV band and in the [0.3 - 0.5] keV, and [0.5 - 1 .0] keV bands 
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Figure 1. Evolution of the EM(T) distributions of the cloud for M — 30 (left panels) and M — 50 (right panels) cases; upper 
panels show the result for models without thermal conduction and radiation (HY models), lower panels for models with both 
effects (RC models). The EM(!T) distributions are sampled every 0.4 t cc since t = 0.2 t cc . The time sequence follows the color 
codes as reported in the color bar. The temperatures of the shocked ambient medium (r ps h, see Tab. [5) and of the shocked cloud 
medium assuming negligible thermal conduction (T sc \, see text) are marked with vertical dashed and dotted lines, respectively. 



typically selected for the analy sis of evolved SNR shock-cloud 
interactions (see, for instance, Miceli et al. The figure 

shows the X-ray luminosity, Lx, of the shocked cloud only and 
does not consider the contribution originating from the shocked 
ambient medium surrounding the cloud. As expected, Lx is 
larger in the hotter Al = 50 case than in the M — 30 case 
in all the energy bands and in both HY and RC models. In all 
the cases, the X-ray luminosity of the shocked cloud reaches its 
maximum quite early, around t ~ r cc , and then decreases (even 
by one order of magnitude, for instance, in HYm50cl0 in the 
[0.3 - 0.5] keV band); therefore, in the X-ray band, shocked 
interstellar gas clouds will be preferentially visible during the 
early phases of the shock-cloud collision. 

We now discuss in detail the evolution of the cloud mor- 
phology as detected in the X-ray band. We expect the rich- 
est scenario from the hottest M — 50 case. Fig.[3]shows 2-D 
sections in the (x, z) plane of the mass density distribution, p, 
the temperature, T, and the X-ray flux in the [0.1 - 10] keV 
band, Fx, derived from HYm50cl0 and RCm50cl0 models at 
t = 1.2 r cc , just after the maximum luminosity of the cloud 
(see Fig. The maps of p and Fx are in log scale to highlight 
structures with very different density and X-ray fluxes. 

At this stage of evolution, the whole cloud material has al- 
ready been shocked. The size of the cloud ( < 1 pc) is smaller 
than that of the original unshocked cloud (2 r c \ = 2 pc) due to 



the cloud compression. The core of the cloud is a high density 
region (an and Or, see upper panel in Fig. [3J where primary 
and reverse shocks transmitted into the cloud are colliding (see 
also Paper I). A low density region (ch and cr) due to a large 
vortex ring has developed just behind the cloud. In HYm50cl0, 
hydrodynamic instabilities are developing at the cloud bound- 
aries: the combined effect of instabilities and shocks transmit- 
ted into the cloud leads to unstable high-density regions at the 
cloud boundaries (go). In RCm50cl0, the thermal conduction 
suppresses the instabilities and leads to smooth gradients of 
density and temperature from the inner part of the cloud to the 
ambient medium. In both cases, the global forward shock has 
converged on the symmetry axis (z-axis), and undergoes a con- 
ical self-re flection, forming the primary Mach reflected shocks 
(da and c/r) and the stem bulge at the base of the secondary 
v ortex sheets near the sym metry axis (en and e^, see Fig. 6 
in Poludnenko et al. 2002 for a detailed description of the flow 
structures developing during the shock-cloud interaction). The 
reflected bow shock is visible as a curved region extending into 
the shocked ISM right below and along the sides of the cloud 
(regions / H and / R ). 

The comparison between upper and lower panels of Fig. 
shows that, in both models, the region with the highest X- 
ray flux is in the core of the cloud (regions oh and or). In 
HYm50cl0, the X-ray image of the shocked cloud has a very 
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Figure 2. X-ray light curve of the cloud in the [0.1 - 10] keV 
(top panel), [0.3 - 0.5] keV (middle panel), and [0.5 - 1.0] 
keV (bottom panel) bands of the shocked cloud derived for the 
Al = 30 and Al = 50 shock cases with (RC models) or without 
(HY models) thermal conduction and radiative cooling. 



sharp boundary and even the hydrodynamic instabilities at the 
cloud boundary are clearly marked in the X-rays (region b^); 
high X-ray flux also originates from the unstable high-density 
regions at the cloud boundaries (region gu). In RCm50cl0, in- 
stead, the emission from region £>r is more diffuse, varying 
smoothly in the radial direction from the center of the cloud. 
Fig. |3 also shows that the X-ray emission density in the re- 
flected bow shock (regions fy and /r) and in the primary Mach 
reflected shocks (c/h and c/r) is slightly higher (by a factor ~ 2 
in HYm50cl0, and by a factor ~ 4 in RCm50cl0) than that 
of the post-shock ambient medium not involved in the shock- 
cloud interaction. On the other hand, the low density region 
(<?h and cr) and the stem bulge (<?h and <?r) are characterized 
by very low X-ray emission. 




-1.5 -1 JO -0.5 
x[pc] 



0.5 III 1.5 
x [pc] 



Figure 3. 2-D sections in the (x,z) plane of the mass density 
distribution (gm crrT 3 ; top panels) in log scale, temperature 
(MK; middle panels), and X-ray emission in the [0.1 - 10] 
keV band (erg s~'; bottom panels), in log scale, derived from 
runs HYm50cl0 (left panels) and RCm50cl0 (right panels) at 
t — 1.2 t cc . Labeled regions in the upper panels are discussed 
in the text. 



Fig. |4] shows the map of X-ray emission in the [0.1 - 10] 
keV band (in linear scale) integrated over 10 pc along the LoS 
(see Sect. Pi for HYm50cl0 and RCm50cl0 at four selected 
epochs around the time of the maximum cloud X-ray lumi- 
nosity (see Fig. |2ji. The superimposed contours are the mass 
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Figure 4. X-ray images in the [0.1 - 10] keV band in lin- 
ear scale derived from models HYm50clO (left panels) and 
RCm50clO (right panels) at the labeled times. Contours are the 
mass density distribution in the (x, z) plane, corresponding to 
logp (gmcirT 3 ) = -23.9, -23.7, -23.5, -23.3. Regions "1" 
and "2" in the third row originate from regions «h and gu, re- 
spectively, shown in Fig. [3] 



density distribution (in log scale) in the (x, z) plane. The third 
row of plots corresponds to the time of Fig.[3](f = 1.2 t cc ). 
The highest integrated emission originates in the core of the 
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Figure 5. As in Fig. |4] for the two labeled times and for the 
[0.3 - 0.5] keV band. 



shocked cloud in both models; two separate high-emission re- 
gions are visible in HYm50clO around t = t cc : the upper 
region (labeled "1" in the third row of Fig. 0} is the high- 
density region a H discussed above (see Fig. |3j, whereas the 
lower one (region labeled "2") originates from the integration 
along the LoS of the emission of high-density unstable re- 
gions at the cloud boundary (region gu in the upper panel in 
Fig. |3- In HYm50clO, the hydrodynamic instabilities are no 
longer clearly distinguishable after integration along the LoS. 
In RCm50clO, the brightest portion of the shocked cloud cor- 
responds to the high-density region («r), appears significantly 
smaller and shorter-lived than in HYm50clO and, in general, 
the cloud surface brightness rapidly approaches the values of 
the surrounding medium. The reflected bow shock (regions /h 
and /r) and the primary Mach reflected shocks (regions du and 
d-g) have a surface brightness more than a decade lower than 
that of the cloud core. 

Fig-H]shows that the highest emission comes from the core 
of the shocked cloud (regions «h and or) in the [0.3 - 0.5] 
keV band, and is maximum at t « t cc for RCm50clO and at 
~ 1 .2 t cc for HYm50c 1 0. In model RCm50c 1 0, the cloud fades 
out earlier than in HYm50clO because of the dissipation by 
thermal conduction. 

In the higher energy [0.5 - 1.0] keV band, shown in Fig. 
[6] the shocked cloud is bright only around r cc . In particular, 
in HYm50clO, a small fraction of the cloud (the high-density 
region an', see Fig. [5} has significant emission only for a very 
short time around 1.2 t cc (see Fig.|6j. In RCm50clO, instead, 
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Figure 6. As in Fig.|5]for the [0.5 - 1.0] keV band. 



the X-ray image of the cloud appears more extended and dif- 
fuse than in HYm50clO, and the cloud is already visible at 
t = 0.2 t cc (not shown in Fig. [6) and remains bright until 
t = 1.4 r cc . This larger-lasting emission is due to the increased 
cloud X-ray emission at higher energies, determined by the 
thermal conduction that heats the cloud material to higher tem- 
peratures. 

For further details on the evolution of the X-ray emission 
in the three bands selected, see the on-line material. 

3.3. Median energy of X-ray photons 

The map of the median energy of X-ray photons (hereafter the 
MPE map) is a practical tool to convey at the same time both 
spatial and spectral information on th e emitting plasma a t high 
resolution jHong et all2004l see also lMiceli et al.l2005l) . 

Fig-0shows the model MPE maps at the times of Figs. [5] 
and [6] obtained from the spectra synthesized in the [0.1 - 10] 
keV band (see Sect. 12. 2> . By comparing the MPE maps in Fig. 
with the X-ray images in Fig. |4] we note that, during the 
whole evolution, the X-ray emission is high where the median 
photon energy £ is low. This result is evident in the £ versus 
F x scatter plot derived for both HYm50cl0 and RCm50cl0 
models at t — 1.2 r cc (Fig. [8j. In fact, most of the brightest 
pixels are in the shocked cloud, i.e. plasma with temperature 
lower than that of the shocked surrounding medium. Fig.[S]also 
shows that, for F x > 0.5 x 10~ 3 erg s" 1 cirT 2 , £ increases with 
Fx indicating that the cloud plasma is far away from pressure 
equilibrium. 
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Figure 7. Median photons energy maps in the [0.1 - 10] keV 
band derived for H Ym50c 1 (left panels) and RCm50c 1 (right 
panels) at the labeled times. Contour plots as in Fig.|4] 
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Figure 8. Median photon energy, £ , versus X-ray flux in the 
[0.1 - 10] keV band, F x , scatter plot derived for HYm50cl0 
(blue) and RCm50cl0 (red) at t = 1.2 t cc . 



Fig. also shows that, in HYmSOclO, the low £ region 
is rather uniform with £ ss 0.2 keV, and its boundaries are 
sharp. In model RCm50cl0, instead, the thermal conduction 
smoothes the energy gradient: in the low £ region, £ increases 
smoothly from the cloud center to the surrounding medium. 
The minimum £ value is higher (~ 0.3 keV) than that in 
HYm50cl0 (~ 0.2 keV) because of the heat conducted into 
the cloud. 
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4. Discussion and conclusion 

We derived the X-ray emission predicted by hydrodynamic 
modeling of the interaction of a SNR shock wave with an inter- 
stellar gas cloud. Our forward modeling allows us to link model 
results to observable quantities and to investigate the observ- 
ability of features predicted by models. 

Our analysis has shown that the morphology of the X-ray 
emitting structures is significantly different from the morphol- 
ogy of the flow structures originating from the shock-cloud in- 
teraction. For instance, the complex pattern of shocks (e.g. ex- 
ternal reverse bow shock, shocks transmitted into the cloud, 
Mach reflected shocks at the symmetry axis, etc.) as well as 
other flow structures (e.g. hydrodynamic instabilities, the stem 
bulge at the base of the secondary vortex sheets near the sym- 
metry axis, etc.) caused by the shock-cloud collision are visible 
in the density maps, but they are never clearly distinguishable 
in X-ray images (cf. upper panel in Fig.[3]and third row of Fig. 
13 at t — 1.2 t cc ). Indeed, the morphology of the X-ray emitting 
structures appears quite simple in all the cases examined. The 
largest contribution to the X-ray emission originates from the 
core of the cloud where primary and reverse shocks transmit- 
ted into the cloud collide. The bright core is surrounded by a 
diffuse and faint region associated with the outer portion of the 
cloud. The X-ray emission varies smoothly in the radial direc- 
tion from the bright core to the surrounding medium. 

The hydrodynamic instabilities, developing at the cloud 
boundaries in models without thermal conduction, are never 
clearly visible in the X-ray band because faint and washed out 
by integration along the LoS. On the other hand, the interac- 
tion of the instabilities with shocks transmitted into the cloud 
produces a bright region with luminosity comparable to that of 
the cloud core (region labeled "2" in Fig.|4j. At variance with 
models including thermal conduction, therefore, in HY models 
two separate bright regions develop inside the shocked cloud. 
This has an important implication on the diagnostics. In fact, in 
Paper I, we have shown that the thermal conduction is very ef- 
fective in suppressing hydrodynamic instabilities: the evidence 
of these instabilities during the shock-cloud interaction would 
be an indication that the thermal conduction is strongly inhib- 
ited (for instance by an ambient magnetic field). Our analy- 
sis points out that, unfortunately, the X-ray band cannot give 
strong indications about hydrodynamic instabilities in any case. 

Our modeling has also shown that shocked interstellar gas 
clouds reach their maximum X-ray luminosity around t ~ t cc . 
The size of the bright region in X-ray maps varies during the 
shock-cloud interaction: the maximum extension is reached at 
epochs < r cc and is always significantly smaller than the orig- 
inal cloud diameter. The light curve of the shocked cloud and 
the evolution of the bright region indicate that shocked clouds 
are expected to be preferentially observed in the X-rays during 
the early phases of shock-cloud collision. 

As an example, in the M — 50 shock case considered here, 
the shocked cloud has total luminosity in the [0.5 - 1.0] keV 
band Lx > 10 33 erg s _1 during the period 0.4 r cc < t < 1.3 r cc 
(see Fig. EJ. Since the XMM-Newton/EPIC-MOS sensitivity 
limit in the [0.5 - 2.0] keV b and is F m \ c , ~ 10~ 14 erg cirT 2 s _1 
for an exposure time of 10 4 s JWatson et all200ll) . our analysis 



suggests that such emission should be detectable as far as ss 30 
kpc. 

The modeling shows that thermal conduction and radiative 
cooling can lead to two different gas phases emitting in differ- 
ent energy bands: a cooling dominated core which ultimately 
fragments into cold, dense, and compact filaments emitting at 
low energies (e.g. optical band), and a hot thermally conduct- 
ing corona emitting at high energies (e.g. soft X-rays). Both 
phases are clearly present in the M — 30 case but only the hot 
one in the M — 50 case because thermal conduction is highly 
effective. As an implication, we expect that the X-ray emission 
morphology and spectrum of the bright cloud region should be 
sensitive to thermal conduction effects. In fact, thermal con- 
duction makes the X-ray bright region smaller, more diffuse, 
and shorter-lived than that expected when thermal conduction 
is neglected. Also, we found that the median photons energy 
of the bright region is higher in models with thermal conduc- 
tion. As a final diagnostic consideration, we note that observing 
smooth gradients of emission and median photon energy would 
indicate that the thermal conduction is efficient. 

The results presented here illustrate the X-ray radiation 
emitted during the shock-cloud collision. Our analysis provides 
a way: 1) to link the features expected to emit X-rays with 
plasma structures originating during the shock-cloud collision, 
and 2) to investigate the effects on the X-ray emission of the dif- 
ferent physical processes at work. These results will be a guide 
for the interpretation of X-ray observations of middle-aged X- 
ray SNR shells whose morphology is affected by ISM inho- 
mogeneities (e.g. the Cygnus Loop, the Vela SNR, G272.2-3.2, 
etc.). However, a more direct comparison of model results with 
supernova remnant X-ray observations requires to include in- 
strumental response and sensitivity and ISM absorption. In a 
companion paper (Orlando et al., in preparation), we will step 
forward to investigate in detail the direct diagnostics and com- 
parison with the data collected with the latest X-ray instru- 
ments (i.e. Chandra, XMM-Newton). 
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